Current-induced magnetoresistance oscillations in two-dimensional electron systems 
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Electric current-induced magnetoresistance oscillations recently discovered in two-dimensional 
electron systems are analyzed using a microscopic scheme for nonlinear magnetotransport direct 
controlled by the current. The magnetoresistance oscillations are shown to result from drift-motion 
assisted electron scatterings between Landau levels. The theoretical predictions not only reproduce 
all the main features observed in the experiments but also disclose other details of the phenomenon. 
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The effect of a strong dc on magnetoresistance 
has long been an outstanding problem in trans- 
port in two-dimensional (2D) electron systems (ESs). 
In the case of microwave-induced magnetoresistance 
oscillations ) 1 ' 2 ' 3 ' 4 ' 5 ' 6 ' 7 ' a finite dc has been shown to sup- 
press the oscillation and eliminate the negative resis- 
tance existing in the weak current limit. 8,9 Quite surpris- 
ingly, in the case without microwave, a relatively weak 
dc can induce substantial magnetoresistance oscillations 
in 2DESs. 

The current-induced magnetoresistance oscillations 
(CIMOs) were observed in differential magnetoresistance 
in high- mobility Hall-bar specimens The oscillation is 
periodic in inverse magnetic field and its period is tun- 
able by the current density. This discovery was later con- 
firmed in highly doped samples and the differential resis- 
tance oscillating with changing current density at fixed 
magnetic field was also detected*^ Measurements were 
recently reported at higher temperatures 1 - 2 - and careful 
studies in ultrahigh mobility samples were carried outjiS 

In this Letter we show that all these observed CIMOs 
can be well explained with the microscopic balance- 
equation scheme for hot-electron magnetotransport di- 
rect controlled by the current ^ 4 ' 15 ' 16 ' 17 

We consider a 2D system consisting of iV s electrons 
in a unit area of the x-y plane. These electrons, inter- 
acting with each other, are scattered by random impuri- 
ties and by phonons in the lattice. There are a uniform 
magnetic field B = (0, 0, B) along the z direction and 
a uniform electric field E in the x-y plane. The non- 
linear steady state magnetotransport of this system can 
be described in terms of the center of mass and rela- 
tive electron variable o 14 ' 15 ' 16 by the following force and 
energy-balance equations^ 



N s eE + iV s e(v x B) + f (v) = 0, 
v • f (v) + w(v) = 0. 



(1) 

(2) 



Here f(v) = fj(v) + f p (v) is the frictional forces due to 
impurity and phonon scatterings, with the impurity part 



fi(v)=]Tq|||t/(q||)| 2 n 2 (q||,q|| 



(3) 



is the electron energy-loss rate to the lattice. In these 
expressions, f7(qii) and M(q) stand for effective impu- 
rity and phonon scattering potentials, II 2 (q||,f2) is the 
imaginary part of the electron density-correlation func- 
tion at electron temperature T in the magnetic field, 
and function A 2 (q, Q) = 2 il 2 (q,| , Q) [n(Q q /T) -n(Q/T e )} 
(n(x) = l/(c x — 1)), also relevant to phonon emission 
and absorption. The effect of interparticle Coulomb in- 
teractions is included in the II2 function to the degree 
of electron level broadening and screening. The remain- 
ing Il2(q||,0) function is that of a 2D electron gas in a 
magnetic field, which can be written in the Landau rep- 
resentation asi 1 - 4 - 

n 2 (q||,fi) = ^^C„, n ^l^/2)n 2 (n,n',n), (5) 



n 2 (n,n',fi) 



de [/(e) - f(s + Q)] 



x ImG„(e + n)ImG„/(e), 



(6) 



w(v) 



q 



Oj|M(q)| 2 A 2 (q,O q + q|| • v) (4) 



where l B = y/l/\eB\ is the magnetic length, C nyn +i{Y) = 
n\[(n + l)\\- x Y l e- Y [L l n {Y)] 2 with L l n {Y) the associate 
Laguerre polynomial, /(e) = {exp[(e — £t)/T e ] + 
is the Fermi function at electron temperature T c , and 
ImG n (e) is the density of states (DOS) of the broadened 
LL n. 

We model the DOS function with a Gaussian-type form 
for both overlapped and separated LLs (e„ = nu> c is the 
center energy of the nth LL and lj c = eB/m is the cy- 
clotron frequency) f 1 "^ 

ImG„(e) = -^Trjir- 1 exp[-2(e - e„) 2 /r 2 ] (7) 

with a magnetic field B 1 / 2 dependent half width T = 
{^aeuJc/irmiJLo) 1 / 2 expressed in terms of /i , the linear 
mobility at lattice temperature T in the absence of the 
magnetic field, together with a broadening parameter 
a to take account of the difference between the trans- 
port scattering time and the broadening-related quan- 
tum lifetime! 4 -! 7 - We will also use a S-independent Y for 
comparison. 

Equations (Q]) and |(3J) are quite general, applicable to 
current-control magnetotransport in any configuration. 
For an isotropic system where the frictional force is in 
the opposite direction of the drift velocity, we can write 
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f(v) = f(v)v/v and w(y) = w(v). In the Hall configu- 
ration with the velocity v in the x direction v = (v, 0, 0) 
or the current densities J x = N s ev and J y = 0, Eq. ^ 
yields the transverse and longitudinal resistivities 

R xy = E y /J x = B/N s e, (8) 
R xx = E x /J x = ~f(v)/(N2 e 2 v), (9) 

and the longitudinal differential resistivity 

r xx = -(df(v)/dv)/(Ny). (10) 

Equation {§]) confirms the relation E y = BJ x /N s e in the 
Hall configuration of nonlinear magnetotransport even 
with an intense dc flowing in the x direction. 

Expressions (J9j) and (flQ|) are particularly convenient 
to deal with current-induced phenomena. Apparently, 
the velocity v, or the current J x , can affect R xx and r xx 
through the qy ■ v factor in the II2 5 ' v ) function. 
Eqs. (O and ^ indicate that in the case of low elec- 
tron temperature (T c <C ep, the Fermi level) and many 
Landau-level occupation, U 2 , f2) is essentially a peri- 
odical function, i.e. n 2 (q|j , 51 + lj c ) = n 2 (q||,f2). There- 
fore, the impurity-induced resistivity R xx would exhibit 
periodical oscillations when changing drift velocity v or 
changing magnetic field 1/B. We introduce a frequency- 
dimension quantity ujj = 2k$v to trace the change of the 
drift velocity v or the current J x = N s ev, and use the 
dimcnsionlcss ratio 

uij 2mkpv l8n m J x . . 

~ c ~ eB ~\~N s ~e?l3 ' ' 

as the control parameter to demonstrate this oscillation, 
which exhibits an approximate periodicity A(u>j/u> c ) ~ 1. 

In addition to the effect discussed above, a finite cur- 
rent J x may heat the electrons and can also affect the 
longitudinal resistivity through the electron-temperature 
change in the n 2 function. Giving the drift velocity v 
or the current J x = N s ev, the electron temperature T e 
is easily determined by the energy-balance equation 
w(v) + vf(v) = 0, and the longitudinal resistivity is then 
obtained directly from Eq. ^j. 

For GaAs-based high-mobility 2DESs at low tempera- 
tures, the dominant direct contribution to the resistivity 
comes from impurity scatterings and / p is negligible. To 
obtain the electron energy dissipation rate w needed for 
determining the electron heating, we consider scatterings 
from bulk longitudinal and transverse acoustic phonons, 
as well as from polar optical phonons with coupling pa- 
rameters taken as typical values of n-type GaAs^ having 
an electron effective mass m = 0.067 m e (m e is the free 
electron mass). 

Figure 1(a) presents the calculated resistivity R xx 
and differential resistivity r xx versus the inverse mag- 
netic field 1/B in terms of ujj/uj c at lattice temperature 
T = 1 K, for a 2D system with electron density N s = 3.0 x 
10 15 m -2 and linear mobility /xq = 2000m 2 /Vs subject to 
three different bias dc current densities J x — 0.40,0.68 
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FIG. 1: (Color online) (a) Resistivity R xx and differential re- 
sistivity r xx versus u>j/uj c at fixed DC current density J x = 
0.40, 0.68 or 1.14 A/m, and (b) R xx , r xx and electron temper- 
ature T e versus ujj/uo c at fixed magnetic field B = 0.08,0.12 
or 0.24 T, for a GaAs-based 2D system having electron den- 
sity iV s = 3.0 x 10 15 m" 2 , linear mobility fi = 2000m 2 /Vs, 
a = 10 at T = 1 K, assuming SR scattering. 
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FIG. 2: (Color online) R xx , r xx and T c versus Wj/iO c at fixed 
magnetic field B = 0.12 T for the same system as described 
in Fig. 1 but subject to SR, BG or RM scattering potentials. 



and 1.14 A/m, which correspond to u>j/2n = 32.3,62.2 
and 103.6 GHz respectively. The elastic scattering is as- 
sumed to be short ranged, and the broadening parameter 
a = 10, i.e. B 1 / 2 -dependent T, ~ 0.62 K at B — 0.08 T. 
Oscillations in resistivity R xx , especially in differential re- 
sistivity r xx show up remarkably, having an approximate 
period A(ujj/Lu c ) ~ 1. The oscillation amplitude decays 
with increasing u>j/u c (reducing B field, due to increas- 
ing overlap of LLs) at fixed bias current but increases 
with increasing bias current density within the J x range 
shown. The maxima (minima) of the differential resistiv- 
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FIG. 3: (Color online) Calculated R xx , r xx , R xx /Ro, Ar xx = 
r xx — Ro and T c for a GaAs-based 2D system with A^ = 
3.0 x 10 15 m -2 at lattice temperature T = 1.5 K subject to a 
mixed SR and BG elastic scattering. The LL width is taken 
to be B 1/2 dependent with F = 0.4K at B = 0.08 T. The 
insects illustrate the case of fixed LL width F = 0.4 K. 

ity t xx locate quite close to (but somewhat lower than) 
the integers (half integers) of LOj/uj Cl while the maxima 
(minima) of the total resistivity R xx are shifted around a 
quarter period higher. These features are in good agree- 
ment with recent experimental findings^ Note that the 
electron temperature T c (not shown) exhibits only a weak 
variation with changing B field at each fixed dc. 

In Fig. 1(b) we plot R XX1 r xx and electron tempera- 
ture T e versus the dc density in terms of Wj /lu c at fixed 
magnetic fields B = 0.08,0.12 and 0.24T for the same 
system. Remarkable R xx and r xx oscillations with ap- 
proximate period A(ujj/uj c ) ~ 1 and maxima (minima) 
positions similar to those in Fig. 1(a) can be seen in this 
current-sweeping figure but here the oscillation decay 
with increasing LUj/uj c is due to increase in the electron 
temperature. At lower u>j/uj c range, the oscillation am- 
plitude of high B-field case is apparently larger than low 
_B-field case when the electron temperature T c is still in 
the range less than or around 10 K. However, in the case 
of B = 0.24 T the oscillation amplitude decays rapidly 
with increasing u>j/u> c due to the rapid increase of elec- 
tron temperature, which rises up to 20 K range around 
LOj/uj c = 2. In comparison, the amplitude decay is much 
slower in the case of B = 0.08 T because of the slow T c 
rise. Furthermore, the periods of R xx and r xx oscilla- 



tions are also somewhat shrunk by the rise of T e , as can 
be seen at higher orders in the B = 0.24 T case. 

Note that though the periods of these resistance os- 
cillations are roughly the same in terms of u)j/w c , their 
amplitude and the detailed behavior depend strongly on 
the form of the scattering potential t/(qy) in Eq. ([3]). To 
have an idea of this scattering potential effect we plot, in 
Fig. 2, R xx , r xx , and T c as functions of ujj/uj c at a fixed 
magnetic field B = 0.12 T for the same 2D system but 
with the dominant elastic scatterings, respectively, due 
to short-range (SR) disorder, charged impurities in the 
background (BG), or ionized impurities locating a dis- 
tance s = 3nm away from the 2D sheet (RM)4 R xx 
and r xx exhibit the strongest oscillations in the case of 
SR potential, with a feature that the second minimum 
of r xx goes deeper into negative direction than the first 
minimum. In the case of BG scattering R xx and r xx 
oscillations, though weaker than those of SR-scattering 
case, still appear quite substantial and the first minimum 
of r xx turns out to be the deepest one of all minima. In 
the case of RM (s = 3nm) scattering, these current- 
induced resistance oscillations, though existing, become 
much weaker than those of SR and BG scatterings. 

Fig. 3 presents the calculated R xx , r xx , R xx /Rq, 
Ar xx — r xx — Rq (Rq is the resistivity at zero dc 
bias) and T e , for another GaAs-based 2D system with 
N s = 3.0 x 10 15 m- 2 , = 1500m 2 /Vs at T = 1.5K, 
focusing on the first oscillation period of uij/u) c . The 
elastic scatterings are assumed due to a mixture of SR 
and BG impurities (with 2:1 contribution ratio to the 
linear mobility). 16 The broadening parameter is taken to 
be a = 3 (S 1 / 2 -dependent T, ~ 0.4 K at B = 0.08 T), a 
somewhat smaller LL width to keep the system mainly 
in the separated LL regime. The insets of Figs. 3a and 
3c show the results using a ^-independent V = 0.4 K. All 
the main features found in the experiment in the case 
of separated Landau levels^ e.g., the dramatic initial 
suppression of the magnetoresistivity with increasing dc 
density, the widening and deepening of the first minimum 
range with increasing B field, a few resolvable oscillations 
of r xx at lower B, dramatic reduction of it at higher mag- 
netic field, etc., are well reproduced. Note that the width 
of the half zero-bias peak in terms of u>j accurately re- 
flects the width of LLs, as can be seen clearly in Fig. 3(a) 
and its inset. This suggests an ideal way to determine 
the width of LLs. 
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